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Abstract 

We study the QCD phase structure in magnetic field H at low temperature 
T. The hadronic phase free energy in a constant homogeneous magnetic field 
is calculated in one-loop approximation of the chiral perturbation theory. The 
dependence of the quark and gluon condensates upon the temperature and 
field strength is found. It is shown that the chiral phase transition order 
parameter (qq) remains constant provided field strength and temperature are 
related via H = const ■ T 2 . 
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1. The investigation of the vacuum state behavior under the influence of the 
various external factors is known to be one of the central problems in quantum 
field theory. In the realm of strong interactions (QCD) the main factors are the 
temperature and the baryon density. At temperatures below the chiral phase tran- 
sition, T < T c , the dynamics of the system is characterized by confinement and 
spontaneous breaking of chiral symmetry (SBCS). At low temperatures, T < T c , 
the partition function of the system is dominated by the contribution of the lightest 
particles in the physical spectrum. In QCD this role is played by the ir - meson 
which is the Goldstone excitation mode in chiral condensate. Therefore the low 
temperature physics (the hadron phase) enables an adequate description in terms of 
the effective chiral theory Q, A very important problem is the behavior of the 
order parameter (the quark condensate (qq)) with the increase of the temperature. 
In the ideal gas approximation the contribution of the thermal pions into the quark 
condensate is proportional to T 2 0, In chiral perturbation theory (ChPT) the 
two - and three-loops contributions (~ T 4 and ~ T 6 correspondingly) into (qq) 
have been found in || and ||. 

The situation with the gluon condensate (G 2 ) = ((G ! ^ i/ ) 2 )is very different. The 
gluon condensate is not an order parameter in phase transition and it does not 
lead to any spontaneous symmetry breaking (SSB). At the quantum level the trace 
anomaly leads to the breaking of the scale invariance and this in turn results in 
nonzero value of (G 2 ). However, this is not a SSB phenomenon and hence does not 
lead to the appearance of the Goldstone particle. The mass of the lowest excitation 
(dilaton) is directly connected to the gluon condensate, m D oc ((G 2 )) 1 / 4 . Thus the 
thermal excitations of glueballs are exponentially suppressed by the Boltzmann fac- 
tor ~ exp{— rrigi/T} and their contribution to the shift of the gluon condensate is 
small (A(G 2 )/(G 2 ) -0.1 % at T =100 MeV) 0. Next we note that in the one-loop 
approximation ChPT pions are described as a gas of massless noninteracting par- 
ticles. Such a system is obviously scale-invariant and therefore does not contribute 
into the trace of the energy- momentum tensor and correspondingly into (G 2 ). As it 
has been demonstrated in || the gluon condensate temperature dependence arises 
only at the ChPT three-loop level. 

Another interesting problem is the study of the vacuum phase structure under 
the influence of the external magnetic field H. Quarks play an active role in shaping 
the QCD vacuum structure. Being dual carriers of both 'color' and 'electric' charges 
they also respond to externally applied electromagnetic fields. The vacuum of strong 
interactions influences some QED processes has been discussed in Ref. 0. The 
behavior of the quark condensate in the presence of a magnetic field was studied in 
Nambu-Iona-Lasinio model earlier [ID|. For QCD, the analogous investigation in the 



one-loop approximations was done in ||11|| . It was found that the quark condensate 



grows with the increase of the magnetic field H in both cases. It implies that a 
naive analogy with superconductivity, where the order parameter vanishes at same 
critical field, is not valid here. The behavior of the gluon condensate (G 2 ) in the 
Abelian magnetic field is also a nontrivial effect. Gluons do not carry electric charge; 
nevertheless, virtual quarks produced by them interact with electromagnetic field 
and lead to the changes in the gluon condensate. This phenomenon was studied in 
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l^j , |13| based on the low-energy theorems in QCD. The vacuum energy density, 
the values of (G 2 ) and (qq) as functions of H have been found in the two-loop 
approximation ChPT in [|13[| . 

The low-energy theorems, playing an important role in the understanding of the 
vacuum state properties in quantum field theory, were discovered almost at the same 
time as quantum field methods appeared in particle physics (see, for example Low 
theorems [T^]).In QCD, these theorems were obtained in the beginning of eighties 



15[| . These theorems, being derived from the very general symmetrical considera- 
tions and not depending on the details of confinement mechanism, sometimes give 
information which is not easy to obtain in another way. Also, they can be used as 
"physically sensible" restrictions in the constructing of effective theories. An im- 



portant step was made in JT6|, where low-energy theorems for gluodynamics were 
generalized to finite temperature case. 

In the present paper the vacuum free energy in magnetic field at finite tempera- 
ture is calculated in the framework of ChPT. The general relations are established 
which allow to obtain the dependence of the quark and gluon condensates on T 
and H. A new phenomenon is displayed, namely the "freezing" of the chiral phase 
transition order parameter by the magnetic field when the temperature increases. 
The physical meaning of this fact is discussed. 

2. The QCD Euclidean partition function in external Abelian field has the 
following form (T = 

■P 




[DB][Dq][Dq]expi -J dx A J y d 3 xC 



Here the QCD Lagrangian in the background field is 

C = + E QhM ~ - B;) + m q ]q, (2) 

"0 q=u,d 

where Q q - is the matrix of the quark charges for the quarks q = (u, d), and for the 
simplicity the ghost terms have been omitted. The free energy density is given by 
the relation (3VF (T, H, m q ) = —InZ. In the chiral limit m q — > Eq. ([!]) yields the 
following expressions for the quark and gluon condensates 

(qq) (T, H) = ^- mq =o , (3) 

(G)(T ' F) - 4 0(1/ gl) • (4) 

The phenomenon of dimensional transmutation results in the appearance of a non- 
perturbative dimensional parameter 

A = Mexp{f° -^U , (5) 
{Ja s (M) f3(a s ) J 

where M is the ultraviolet cutoff, a s = g 2 /ATT, and f3(a s ) = da s (M)/d InM is the 
Gell-Mann-Low function. In chiral limit (m q = 0) the system described by the 
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partition function ([!]) is characterized by three dimensionful parameters M, T and 
H. The free energy density is renorm-invariant quantity and hence its anomalous 
dimension is zero. Thus F has only a normal (canonical) dimension equal to 4. 
Making use of the renorm-invariance of A, one can write in the most general form 

^ = AV (6) 



,A A 2 

where the function / is still unknown. From ([5]) and (|6|) one gets 

dF OF dA 2 47ra w 4 _ T < 9 ou 9 )F 

With the account of (f|) the gluon condensate is given by 

1 (\ttcv 2 f) f) 

At T = 0, H = we recover the well known expression for the nonperturbative 
vacuum energy density in the chiral limit. In the one- loop approximation ((3 = 
-b a 2 s /2TT, b = (llAe - 2iV»/3) it has the form 

£ „ = F(0,0) = - I A_(G 2 ) (9) 

Let us now derive the low-energy theorems at finite temperature in the presence 
of magnetic field. Strictly speaking, /3-function depends on H and the low-energy 
theorems acquire electromagnetic corrections oc e 4 . However, since the free energy 
is independent of the scale M at which the ultraviolet divergences are regulated, 
we can choose M 2 3> H,T 2 ,A 2 . Hence, we can take the lowest order expression for 
/3-function (/3(a s ) = — b Q a 2 /2tt) and the electromagnetic corrections vanish. Let us 
introduce the field <t(t = x±, x) and operator D, 



^,x) = -^(G^(r,x)) 2 , (10) 



d d 

D = A — T— 2H-. (11) 

dT OH K ' 

Differentiating @ n times with respect to (1/<7q) an d taking into account (0), QTOf) 
and (|ITD one obtains 

D n+l F = t> n (a -(0,0)) 

= J dr n d 3 x n ... J dT 1 d 3 x 1 (a(r n ,x n )...a(r 1 ,x 1 )a(0,O)) c . (12) 

The subscript c means that only connected diagrams are to be included. Proceeding 
in the same way, it is possible to obtain the similar theorems for renorm-invariant 
operator 0(x) which is built from quark and/or gluon fields 

( d d y 
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= J dr n d 3 x n ... J dr 1 d 3 x 1 (a( y T n ,x n )...a(ri,x 1 )O(0,0)) c . (13) 

where d is the canonical dimension of operator O. If operator O has also anomalous 
dimension, the appropriate 7-function should be included. 

3. The above equations enable to obtain the values of the condensates as func- 
tions of T and H provided the free energy density is known. To get the latter the 
ChPT will be used. At low temperatures T < T C (T C is the chiral phase transition 
temperature) and for weak fields H < n\ adr ~ (inF^) 2 the characteristic momenta 
in the vacuum loops are small and theory is adequately described by the low-energy 
effective chiral Lagrangian L e jf 0, ||. This Lagrangian can be represented as a 
series expansion over the momenta (derivatives) and quark masses 

L eff = + + L (6) + ... (14) 



The leading term in (14) is similar to the Lagrangian of the non- linear sigma model 
in the external field 

F 2 

L {2) = -jTr(V ^U + V M C7) + ERe Tr(MU + ), 

V,JJ = d li U-i[U,V lx \. (15) 

Here U is a unitary SU(2) matrix, F n = 93MeV is the pion decay constant, and 
E has the meaning of the quark condensate S = = \(dd)\. The external 

Abelian magnetic field H is aligned along the z -axis and corresponds to V^{x) = 
(r 3 /2)A fl (x) with the vector-potential chosen as A^(x) = S^Hxi. The mass 
difference between the u and d quarks appears in the effective chiral Lagrangian 
only quadratically. Further, to obtain an expression for the quark condensate in the 
chiral limit we use only the first derivative with respect to the mass of one of the 
quarks. Therefore, we can neglect the mass difference between the u and d quarks 
and assume the mass matrix to be diagonal M = ml. 
At T < T c , H < 

l^hadr ^ ne QCD partition function coincides with the partition 
function of the effective chiral theory 

Z eff [T, H] = e-^ffFM = Zo [H] J [DU] exp{- f dx 4 J d 3 xL eff [U, A}} (16) 

At the one- loop level it is sufficient to restrict the expansion of L e // by the quadratic 
terms with respect to the pion field. Using the exponential parameterization of the 
matrix U(x) = exp{iT a ir a (x)/F 7T } one finds 

L {2) = \{d^) + l -Ml(ixy + {d^ + + iA^ + )(d,7t- - iA.tt-) + M 2 n+n-, (17) 
where the charged 7r =l= and neutral tt° meson fields are introduced 



7T 



± -(7r 1 ±m 2 )/v / 2, n° = n 3 
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Thus (pi]) can be recasted into the form[| 

Z* f [T, H] = Z p jZ [H] J [Dn°] [Dn+] [Dir~] exp{- J* dx, ^ d 3 xL^ [tt, A]} (19) 

where partition function is normalized for the case of perturbation theory at T = 
0,H = 

Z p ,. = [det(-dl + Ml)}~^. (20) 
Integration of fll9|) over 7r-fields leads to 



Zftt = Z p jZ [H][det T (-dl + M^)]-^[det T (-|^i 2 + M 2 )] -1 , (21) 

where = <9 M — iA^ is a covariant derivative and a symbol T means that the 
determinant is calculated at finite temperature T according to standard Matsubara 
rules. Taking (p0[) into account and regrouping multipliers in fl2T|) one gets the 

7R 
J eff 



following expression for Z^ 



Z* f [T,H] = Z [H] 



detrHg + Ml) 
det(-dl + M%) 



-1/2 



det(-|^| 2 + M|) 

det( _ 9 2 +M 2 ) 



detr(— 




2 + Ml) 


det(- 




2 + Ml 



-I -1 



Then the effective free energy can be written in the form 

F e R ff (T, H) = ~]n Z* SJ = ^ + FAT) + (H) + F S (T, H). 
Here F T o is the free energy of massive scalar boson 

F n0 (T) = T j-^M 1 - ex P(- VP 2 + M '/ T ))> 



(22) 



(23) 



(24) 



F n ± is a Schwinger result for the vacuum energy density of charged scalar particles 
in the magnetic field. 



FAH) 



167T 2 JO S 



00 ds 



-Mis 



Hs 



smh.(Hs) 



1]. 



(25) 



and 



Fs(T,H) = — Y,J o dkln(l-exp(-u n /T)), 



71- 



n=0 1 



0O r . 



k 2 + Ml + if (2n + 1), 



(26) 



1 The partition function -Z^j describes charged ^ and neutral n° ideal Bose gas in magnetic 
field. Relativistic charged Bose gas in magnetic field at finite temperature and density with ap- 
plication to Bose-Einstein condensation and Meissner effect was studied in Refs. [jl7|, pq |, ]i~9[ ], 
121. 
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where uj n are Landau levels of the 7r =1= mesons in constant field H . | 

4. The free energy F^jj- determines the thermodynamical properties and the phase 
structure of the QCD vacuum state below the temperature of the chiral phase tran- 
sition, i.e. in the phase of confinement. 

Equations @ and (^) describe the dependence of (G 2 ) on T and H. The action 
of the operator D on F$ f leads to DF^o (T) = since M 2 — > and F w o(T) ~ T 4 
in chiral limit and thus (4-Td / 1 dT)F^o{T) = 0. It can be easily shown by direct 
calculation that DF S (T,H) = 0. The nontrivial dependence of (G 2 ) on H arises 
only due to Schwinger term F n ±(H) 

a 2 

(G 2 )(T,H) = (G 2 ) + -^—H 2 (27) 

Next we note that because of the asymptotic freedom the QCD /3-function, (3{a s ) = 
—boa 2 /2ir + ... is negative and hence the gluon condensate diminishes with the H 
increasing 

(G 2 )(T,H) = (G 2 )-^-H 2 . 

3&o 

Thus, the temperature corrections to the gluon condensate in magnetic field vanish 
at the ChPT one-loop level. 

In order to get the dependence of the quark condensate upon T and H use is 
made of the Gell-Mann-Oakes-Renner relation (GMOR) 

F 2 M 2 = --(m u + m d )(uu + dd) = 2mS (28) 

Substituting (|3|) into (|3]), calculating the derivative over M 2 and then taking the 
limit M 2 — > one gets 



mT,H) = (*>(! - \ • J + ^ln2 - ^(^)) 

<PW = E / 7\ < yv 1V ^ = ^ + 2n + 1 (29) 

r^Vo tu„ x exp AcUr, uc — 1 



n=0 



Now we consider various limiting cases. In the strong field, y H ^> T (A ^> 1), the 
lowest Landau level (n = 0) gives the main contribution to the sum ( f29|) 



V?(A»l) = ^e- A + 0( e ^ A ). (30) 



In the opposite limit of weak field, y H <C T(A <C 1), the sum in (p9|) is calculated 



with required accuracy using the Euler-MacLaren formula. Furthermore, one gets 

technically, a transition for the free energy F = iTr ln(p|+o;§(p)) from the vacuum case (H = 
0, T — 0) to the case of H ^ 0, T is straightforward. Omitting the details of the calculations, we 
note that, eventually, this transition reduces to the substitutions p± — > LOk = 2irkT (k = 0, ±1, ...), 
u a = VF+Mj - w? = y<pl+MZ + H(2n+l) and Tr - ^ £~ E^-oo /-» ff > where 
the degeneracy multiplicity of Hjlix has been taken into account for the Landau levels. Performing 
summation over Matsubara frequencies, we obtain (26). 
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V(A « 1) = -T2 + 7T7T + 7 ln A + C + ^ A 2 + 0(A 4 ), (31) 



the following result with the use of the asymptotic expansion of integral (|29|) [ gj] at 
A<1 

^ 2 1 7tt 1 1, , ^ C(3) 
— -Tg + 7TTT + 7 In A + C + ^4 
6 A 2 24 A 4 487r 2 

here C = J (7 - ln4vr - §), 7 = 0.577... is Euler's constant and C(3) = 1.202 is 
Riemann zeta function. Thus, one obtains the following limiting expressions for the 
quark condensate in the chiral limit in a magnetic field at T 7^ 

(q q )(T,H) IT 2 H E ^ 

(qq) 3 8F 2 + (4vrF.) 2ln2 (2tt)V**? C ' ^ >> (32) 

and 

(qq)(T,H) _ T 2 H iVHT H H r— 

(qq) ~ 1 ~ 8F 2 + JWJ 2A ~ 48^ " JwJ 2 ^ T 2 ' ^ (33) 

where A = In 2 + 8C ~ 4.93. 

In the framework of ChPT the quark condensate fl29|) at if 7^ 0, T 7^ is de- 
termined by three dimensionless parameters H / '(AnF^) 2 , T 2 /F 2 and A = \fH/T. 
The quantity A is a natural dimensionless parameter in this approach. The motion 
of a particle (massless pion) in the field is characterized by the curvature radius of 
it's trajectory, and in the magnetic field this is the Larmor radius Rl = 1/y/H. 
On the other hand, there is another length It = l/T - "temperature length" at 
T 7^ 0. Therefore, charged 7r ± -mesons in magnetic field effectively acquire "mass", 
m e ff = \/~H, determined by the lowest Landau level, when Larmor radius of a 
particle in the field is much less than i^(A ^> 1). Correspondingly, their contri- 
bution to the shift of the chiral condensate is suppressed by the Boltzman factor 
oc exp{— m e ff/T}. In the weak field limit 7r ± -mesons give standard temperature 
one-loop approximation ChPT contribution to (qq). Besides, additional tempera- 
ture and magnetic corrections appear. Neutral 7r°-meson contributes to (qq) (T, H) 
as usual massless scalar particle. 

An interesting phenomenon reveals itself in the vacuum QCD phase structure 
under consideration. One can find from ( P9"[ ) such a function H(T) that the chiral 
condensate (qq){T, H) remains unchanged when the temperature and magnetic field 
change in accordance with H* = H(T). Then if* is found by solving the following 
equation (see (g|) (qq)(T,H*) - (qq) = 0) 

l-- 3 T A 2 ln2 + ^AV(A) = (34) 



The numerical solution of (|34D yields A* = 0.111... Thus, quark condensate stays 
unchanged when T and if are increased according to H = 0.013 • T 2 . Hence it 
is possible to say that the order parameter (qq) of the chiral phase transition is 
"frozen" by the magnetic field. 

Note that H (T c ) / (An F n ) 2 ~ 2 ■ 10" 4 < 1 at T = T c ~ 150MeV and therefore 
the above relations remain valid up to the deconfined phase transition point. In the 
vicinity of T c the effective low energy chiral Lagrangian fails to provide an adequate 
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Figure 1: Quark condensate S(T, H)/H as function of x — T '/ 'y/8F n and y = 

description of the QCD vacuum thermodynamical properties, and strictly speaking 
becomes physically invalid. The following is worth noting. In deriving (|29|) , at the 
first step the physical quantity as functions of M n where obtained, and only then 
the chiral limit M v — > was taken. Acting in the inversed sequence we would have 
obtained all temperature corrections to condensate identically equal to zero. This 
points to the fundamental difference of the two cases: the exactly massless particle 
and the particle with infinitesimal small mass. 

5. It has been shown in the present letter that the quark condensate is "frozen" 
by the magnetic field when both temperature T and magnetic field H are increased 
according to the H = const ■ T 2 law. This points to the fact that the direct anal- 
ogy between the quark condensate in QCD and the theory of superconductivity is 
untenable. In the BCS theory the Cooper pairs condensate is extinguished by the 
temperature and magnetic field. The "freezing" phenomenon can be understood in 
terms of the general Le Chatelier-Braun principle 0.The external field contributes 
into the system an additional energy density H 2 /2. The system tends to compensate 
this energy change and to decrease the free energy by increasing the absolute value 
of the quark condensate: Ae v = — m\L(H) — S(0)| < 0. On the other hand, if the 
temperature of the system is increased (by bringing some heat into it) the processes 
with heat absorbtion by damping the condensate are switched on. The interplay 
of these processes is at the origin of the above "freezing" of E(T, H). Next, since 
gluons do not carry electric charge, the magnetic field affects the gluon sector of the 
vacuum only indirectly via the quark sector and thus the Le Chatelier-Braun prin- 
ciple is not applicable directly to the gluon condensate. For the same reason gluon 
condensate decreases nonlinearly with H increasing according to A(G 2 ) oc — H 2 , 
while for the quark condensate AS oc H. 

The author is grateful to B.L.Ioffe, V.A.Novikov, Yu.A.Simonov, A.V.Smilga and 
S.M.Fedorov for comments and discussions. The financial support of RFFI grant 
00-02-17836 is gratefully acknowledged. 

3 The external action disturbing the system from the equilibrium state induces processes in this 
system which tend to reduce the result of this action 
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